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ABSTRACT 

Histone variants play an important role in shap- 
ing the mammalian epigenome and their aberrant 
expression is frequently observed in several types 
of cancer. However, the mechanisms that mediate 
their function and the composition of the variant- 
containing chromatin are still largely unknown. A 
proteomic interrogation of chromatin containing the 
different H2A variants macroH2A.1.2, H2A.Bbd and 
H2A revealed a strikingly different protein compo- 
sition. Gene ontology analysis reveals a strong en- 
richment of splicing factors as well as components 
of the mammalian replisome in H2A.Bbd-containing 
chromatin. We find H2A.Bbd localizing transiently 
to sites of DNA synthesis during S-phase and dur- 
ing DNA repair. Cells that express H2A.Bbd have a 
shortened S-phase and are more susceptible to DNA 
damage, two phenotypes that are also observed in 
human Hodgkin's lymphoma cells that aberrantly ex- 
press this variant. Based on our experiments we con- 
clude that H2A.Bbd is targeted to newly synthesized 
DNA during replication and DNA repair. The transient 
incorporation of H2A.Bbd may be due to the intrin- 
sic instability of nucleosomes carrying this variant 
or a faster chromatin loading. This potentially leads 
to a disturbance of the existing chromatin structure, 
which may have effects on cell cycle regulation and 
DNA damage sensitivity. 



INTRODUCTION 

Chromatin consists of an equal mass of DNA and histone 
proteins and is bound by a plethora of associated factors (1- 
6). The fundamental unit of chromatin is the nucleosome, a 
structure that contains 147 bp of DNA wrapped around a 
histone octamer made up by the four core histones H2A, 
H2B, H3 and H4. Individual nucleosomes can differ from 
each other by a multitude of post-translational modifica- 
tions and by the incorporation of distinct histone variants. 
Whereas relatively few variants for H4 and H2B are known, 
most eukaryotes have several isoforms of H3 and H2A (7). 
Canonical histones, which constitute the majority of the hi- 
stones in a cell, are incorporated during S-phase whereas 
most histone variants are assembled independently of DNA 
synthesis. Genome-wide association studies in various or- 
ganisms demonstrate a non-random association of histone 
variants across the genome (8) suggesting a possible func- 
tion in setting up defined chromatin structures. 

The structure of chromatin can undergo massive and 
global rearrangements during cellular differentiation (9- 
12), which are often accompanied by a substantial differ- 
ence in the chromatin composition. Changes in chromatin 
structure and nuclear morphology are frequently observed 
during oncogenic transformation and are often used as di- 
agnostic markers. In this light it is noteworthy that the over- 
expression of specific histone variants has recently been 
shown to be a hallmark of several cancer types. The gene 
encoding the replication-independent H3.3 variant for ex- 
ample has been found to carry a driving mutation in paedi- 
atric glioblastoma (13,14). The two prevalent mutations re- 
sult in a disturbed modification of H3K27 and/or H3K36 
presumably causing a mis-regulation of epigenetic marks. 
The centromeric H3 variant CENP-A on the other hand 
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is overexpressed in primary colorectal cancers, which leads 
to a mistargeting of this variant (15) resulting in mitotic 
defects and an increase in aneuploidy, a hallmark of can- 
cerous cells (16). Besides H3, several H2A variants have 
also been shown to be mis-regulated in cancer (17-20). The 
two most divergent H2A variants H2A.Bbd and macroH2A 
show a striking dichotomy when it comes to their involve- 
ment in cancer: whereas H2A.Bbd RNA is overexpressed 
in strongly proliferating Hodgkin's lymphoma (HL) cells 
(17), the expression of macroH2A significantly suppresses 
melanoma progression (18). These findings suggest that the 
two variants may have opposing properties in regulating cel- 
lular physiology. 

MacroH2A has been named after its large C-terminal 
macro domain (21). Initially it was shown to be enriched 
on the inactive female X-chromosome in mammals (22,23) 
but has later been shown to also localize to large autoso- 
mal domains (1-3,8,24-26). MacroH2A stabilizes nucleo- 
some core particles (27) and interferes with the remodeling 
activity of the SWI/SNF complex (22,28). Consistent with 
a rather repressive function, macroH2A mostly localizes to 
heterochromatic regions and correlates with H3K27me3- 
rich regions (24,29,30). H2A.Bbd preferentially binds H4 
hyperacetylated regions and was named Barr body-deficient 
according to its striking exclusion from the inactive fe- 
male X-chromosome when expressed ectopically (17,31). 
H2A.Bbd-containing nucleosomes organize less DNA and 
are less stable than canonical nucleosomes (16,32). The 
mouse genome contains four genes encoding putative ho- 
mologues of H2A.Bbd termed H2A.Lapl-4 (33,34). The 
H2A.Lapl isoform is expressed in testis and brain and lo- 
calizes to the transcription start sites of active genes in 
round spermatids where it can even localize to the inactive 
X-chromosome (27,29,30,35). When ectopically expressed 
in HeLa cells, H2A.Bbd is enriched on active genes where 
it is involved in transcriptional activation and mRNA pro- 
cessing (29,30,36). This potential role of H2A.Bbd in tran- 
scriptional activation and splicing regulation may have im- 
portant implications in HL cells that express increased lev- 
els of H2A.Bbd and also show an increased mobilization of 
transposable elements (37,38). 

Despite a large body of work with regard to the struc- 
ture and localization of nucleosomes carrying these two hi- 
stone variants, surprisingly little is known about their func- 
tion under physiological or pathophysiological conditions. 
Mice lacking both macroH2A genes are viable and fertile 
(39-41) but their somatic cells are more susceptible to epi- 
genetic remodeling compared to wild type suggesting that 
this histone variant acts as a barrier to induced pluripotency 
(42). Mice lacking the H2A.Bbd homologues H2A.lapl- 
4 have not been reported so far but the testis-specific ex- 
pression pattern suggests a role of this histone variant dur- 
ing male gametogenesis (35,43). In order to better under- 
stand the role of H2A.Bbd upon expression in somatic cells 
and its role during tumor progression, we ectopically ex- 
pressed the H2A.Bbd variant in mouse embryonic fibrob- 
lasts (MEFs) and analyzed the composition of chromatin 
domains containing H2A.Bbd by comparing it to domains 
that either contained canonical H2A or the differentially lo- 
calizing macroH2A.1.2 variant. 



MATERIALS AND METHODS 

Cell culture and cell transfection 

Female MEF cells were grown under standard condi- 
tions in Dulbecco's Modified Eagle Medium (DMEM) 
supplemented with 10% fetal bovine serum (Sigma) and 
penicillin/ streptomycin (PA A). Human HL cell lines L- 
1236 (44), L-540 (45), L-428 (27,46), HDLM-2 (28,46) 
and KM-H2 (24,30,47) were obtained from the Deutsche 
Sammlung fur Mikroorganismen und Zellkulturen (Braun- 
schweig, Germany) HL cell lines and control Epstein- 
Barr virus-immortalised cell lines (LCL) were cultured 
in RPMI 1640/GlutaMAX (Gibco) supplemented with 
10% fetal bovine serum (Sigma), betamercaptoethanol and 
penicillin/ streptomycin (PAA). Peripheral blood mononu- 
clear cells (PBMCs) from healthy donors and patients with 
HL were isolated as described (31,48) with informed con- 
sent and approval by the ethics committee of the Medical 
Faculty of the Martin Luther University Halle- Wittenberg. 
MEF cells were transfected with GFP-H2A.Bbd and GFP- 
H2A expression plasmids using FuGENE® HD transfec- 
tion reagent (Roche) according to the manufacturer's in- 
structions. Stable cell lines were selected using G418 (PAA) 
at a concentration of 600 |xg/ml. The GFP-macroH2A ex- 
pression plasmid was transfected using the TransIT LT1 
reagent (Mirus) according to the manufacturer's instruc- 
tions. For live cell imaging, MEF cells stably expressing 
H2A.Bbd were transfected with mRFP-PCNA (49) using 
Amaxa nucleofection as previously described (50). 

FACS sorting and FACS analysis 

GFP single cell sorting was carried out using a MoFlo 
cell sorter (Beckman Coulter). Single cell clone suspensions 
of GFP-H2A.Bbd, GFP-macroH2A and GFP-H2A MEFs 
were analyzed using a FACS Canto flow cytometer (Beck- 
man Coulter). The GFP signal was detected in the Fluo- 
rescein Isothiocyanate (FITC) channel, and the data analy- 
sis to quantify the GFP positive cells was performed using 
Flow Jo software. 

Chromatin purification 

Adherently growing MEFs were detached by trypsinization 
followed by centrifugation at 200 x g. Cells were washed 
twice with phosphate buffered saline (PBS) and nuclei were 
isolated by the addition of ice-cold PBS/0.3% Triton-X. Af- 
ter a 10 min incubation step at 4°C, nuclei were pelleted for 
10 min at 1000 x g, washed once with PBS and resuspended 
in EX-100 (10 mM Hepes pH 7.6, 100 mM NaCl, 1.5 mM 
MgCl 2 , 0.5 mM EGTA, 10% v/v glycerol) freshly supple- 
mented with proteinase inhibitors. Chromatin was solubi- 
lized by 20 min micrococcal nuclease (MNase) digestion 
at 26° C in the presence of 2 mM CaCl2. The reaction was 
stopped by adding l/50th volume of 0.5 M EGTA. All fur- 
ther steps were carried out at 4°C. The first supernatant (SI) 
was recovered after centrifugation at 1000 x g and the pel- 
let was resuspended in resuspension buffer (lx PBS, 150 
mM NaCl, 2 mM EDTA, 0.1% Triton) and over-head ro- 
tated over night (o/n) to further extract solubilized chro- 
matin fragments. The second supernatant (S2) was recov- 
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ered by a 30 min centrifugation, pooled with the first super- 
natant to yield the purified chromatin. A chromatin aliquot 
was used to check the MNase digestion by DNA purifica- 
tion and agarose gel electrophoresis to measure the size of 
chromatin fragments. 

GFP pull-down 

MNase-digested chromatin was precleared with 20 \x\ of a 
1:1 (v/v) slurry of agarose beads in EX- 100 (10 mM Hepes 
pH 7.6, 100 mM NaCl, 1.5 mM MgCl 2 , 0.5 mM EGTA, 
10% v/v glycerol) for 1 h at 4°C. The supernatant was col- 
lected after 2 min of centrifugation at 2000 x g, 4°C. Twenty 
microliters of a 1 : 1 slurry (v/v) of GFP-binder (Chromotek) 
were equilibrated three times with EX- 100 and then added 
to the supernatant. Immunoaffinity purification was per- 
formed for 2 h at 4°C. The beads were washed three times 
with resuspension buffer (lx PBS, 150 mM NaCl, 2 mM 
EDTA, 0.1% Triton) plus proteinase inhibitors and cen- 
trifuged at 4°C at 2000 x g. Fifteen microliters of 2x 
Lammli buffer were added to the beads and boiled for 15 
min at 95°C . Samples were run on a 15% SDS gel for 2 
cm and stained with Coomassie Brilliant Blue G250 for 20 
min. After destaining with 10% acetic acid, every lane was 
cut into eight slices and every slice was collected into a tube. 

Trypsin digestion and mass spectrometry 

Gel pieces were washed once with milliQ water and two 
times with 20 mM NH4HCO3. Destain of the gel pieces was 
preformed in 50% ACN, 50% NH4HCO3 at 37°C for 1 h. 
Samples were washed two times with 20 mM NH4HCO3 
and dehydrated with 100% ACN (supplier and quality grade 
ultrapure). Reduction was performed in DTT 10 mM di- 
luted in 20 mM NH4HCO3 for 1 h at room temperature. 
Alkylation of the gel pieces was done in 50 mM IAA diluted 
in 20 mM NH4HCO3 at room temperature in the dark for 30 
min. Before tryptic digestion, the samples were washed once 
with 20 mM NH4HCO3 and dehydrated with 100% ACN. 
30-50 |ul1 of trypsin (25 ng/ (xl) was added to the gel pieces 
in ice till the gel pieces hydrated completely. Not absorbed 
trypsin was removed and 50 (Jul of 20 mM NH4HCO3 was 
added. The samples were incubated overnight at 37°C. The 
supernatant containing peptides was collected and placed in 
a fresh vial. Gel pieces were incubated twice in 50% ACN, 
0.25% TFA, 10 min at room temperature followed by two 
times 100% ACN, 10 min at room temperature. All super- 
natants were pooled, dried and then resuspended in 15 |xl 

0. 1% FA. Fifty percent of the sample was injected in an Ul- 
timate 3000 HPLC system (LC Packings Dionex). Samples 
were desalted on-line by a CI 8 trapping column (300 |xm 

1. d. 5 mm, packed with C18 PepMapTM, 5 |xm, 100 A by 
Thermo-Fisher Scientific), and subsequently peptides were 
separated with a gradient from 5% to 60% acetonitrile in 
0.1% formic acid over 40 min at 300 nl/min on an analyti- 
cal column (75 |xm i.d. x 10 cm, packed with C18-PepMap, 
3 |xm, 100 A). The effluent from the HPLC was directly 
electro-sprayed into the LTQ Orbitrap mass spectrometer 
(Thermo Fisher Scientific). The mass spectrometry (MS) in- 
strument was operated in the data-dependent mode to au- 
tomatically switch between full scan MS and MS/MS ac- 



quisition. Survey MS spectra (m/z 350-2000) were acquired 
in the Orbitrap with resolution 60 000 at m/z 400. Frag- 
ment ion spectra were recorded in the linear trap of the in- 
strument. For all measurements with the Orbitrap detector, 
three lock-mass ions from ambient air (m/z = 371.10123, 
445.12002, 519.13882) were used for internal calibration as 
described. The six most intense peptide ions with defined 
charge states of 2+ and higher were sequentially isolated in 
a 2 Da window around the precursor signal (window = 350- 
2000 m/z) to a target value of 10 000 and fragmented in 
the linear ion trap by collision-induced dissociation (CID). 
The ion selection threshold for triggering acquisition of an 
MSMS spectrum was 10 000 counts. Previously fragmented 
precursor signals were excluded from repeated fragmenta- 
tion analysis for 30 s. Typical mass spectrometric condi- 
tions were: spray voltage, 1.4 kV; no sheath and auxiliary 
gas flow; heated capillary temperature, 200° C; normalized 
collision energy, 35% for CID in linear ion trap. 

Mass spectrometry data analysis 

For protein identification, the raw data were analyzed 
with the Andromeda algorithm of the MaxQuant protein 
analysis package (version 1.2.2.5) against the IPI mouse 
v.3.68.Fasta including reverse sequences and contaminants. 
The trypsin/P enzyme was selected, allowing for maximally 
two missed cleavages. Carbamidomethylation of C was set 
as fixed modification whereas methionine oxidation and 
protein N-acetylation were included as variable modifica- 
tions. The mass tolerance of the initial search was 20 ppm; 
after recalibration, a 6 ppm mass error was applied for the 
main search. Fragment ions were searched with a mass off- 
set of 0.5 Da using the six most intense signals per 100 
Da. Searching for secondary peptide hits within already as- 
signed MSMS spectra was enabled. The search results were 
filtered with a peptide and protein false discovery rate of 
0.01 with a minimum peptide length of six amino acids. 
Protein identifications with at least one unique or razor 
peptide were accepted. For the quantification, the Intensity 
Based Absolute Quantification (iBAQ) values were calcu- 
lated from peptide intensities and the protein sequence in- 
formation (51) of unmodified, M/ oxidated and acetylated 
peptide species with a minimum of two peptides per pro- 
tein. The mass spectrometry proteomics data have been de- 
posited to the ProteomeXchange Consortium (http://www. 
proteomexchange.org) via the PRIDE partner repository 
(52) with the data set identifier PXD000789. 

Statistical data analysis 

Statistical analysis of protein enrichment was performed 
from the three affinity purification- MS replicates of the 
three GFP-H2A variants. First, all replicates were normal- 
ized using the median of abundant histone H4 peptides, in- 
cluding only peptides without modifications, less than two 
missed cleavages and excluding tryptic cleavage after Pro, 
as normalization factor. No-GFP data were not normal- 
ized, since detected H4 levels were too low. Second, normal- 
ized iBAQ values were log2 transformed and missing values 
were imputed from a Gaussian distribution that was shifted 
by a difference of about 2.5 and had a width of 0.3 using 
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Figure 1. Female MEF single clones stably express the GFP-tagged histone variants. (A) FACS profile of female MEF cell populations derived from 
sorted single cell clones that express H2A.Bbd, macroH2A.1.2 and H2A tagged with GFP (solid line) in comparison to untransfected control female MEF 
cells (dotted and filled histogram). (B) Western blot analysis of GFP-tagged histone variant expression levels with antibody against GFP. Detection of 
histone H3 served as loading control. (C) Nuclear localization of the histone variants GFP-H2A.Bbd and -macroH2A.1.2 in female MEFs stained with 
H3K27me3 (red). Arrows indicate the Barr body. Scale bar: 5 |jim. 

the Perseus software tool (53). The average difference of 
the GFP-H2A variants over cells not expressing GFP (no- 
GFP) was calculated and only proteins having a difference 
of more than 3.5 and quantitative information in at least 
two experiments were kept. To determine enrichment fac- 
tors and /^-values for enrichment, ANOVA was performed 
in DanteR and resulting ^-values were subsequently cor- 
rected for multiple hypothesis testing using the Benjamini- 
Hochberg method in DanteR (54). Data were plotted with 
Perseus software. The gene ontology was performed using 
the Database for Annotation, Visualization and Integrated 



Discovery (DAVID) version 6.7. (55,56) Functional Anno- 
tation Charts were created using default settings. 

Immunofluorescence 

Cells grown on glass coverslips were fixed in PBS/ 3. 7% FA 
(Roth:7398.1) for 10 min at room temperature. After fixa- 
tion, the cells were washed in PBS and then permealized in 
3.3 mM sodium citrate and 0.1% Triton X-100 for 5 min at 
room temperature. The cells were then washed two times 
with PBS and two times with washing solution I (PBS- 
Tween (0.1%) plus 0.25% BSA). Blocking was performed 
with blocking solution I (5 ml washing solution 1/2.5% 
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Figure 2. Chromatin containing GFP-tagged histone variants is released after MNase digestion. (A) Chromatin purification protocol (left) and DNA 
analysis of the first and second supernatants containing chromatin (SI and S2, respectively). (B) Western blot analysis of the GFP-tagged histone variant 
release compared to H3 in presence or absence of chromatin MNase digestion. 



BSA) for 30 min at room temperature. Primary antibodies 
were diluted in blocking solution I. Incubation with the pri- 
mary antibody was performed overnight at 4°C. The next 
day after washing with washing solution I the cells were 
incubated with the secondary antibodies diluted in block- 
ing solution II (blocking solution I plus 10% normal goat 
serum) for 1 h at room temperature. Slides were mounted 
with Vectashield containing 4 / ,6-diamidino-2-phenylindole 
(DAPI; Vectore Laboratories) and sealed with nail polish. 
Images were acquired using a Zeiss LSM 510 META con- 
focal microscope equipped with a 63 x oil immersion ob- 
jective NA = 1 .4 or a Leica TCS SP5 confocal microscope 
with a 63 x glycerol immersion objective NA =1.3. Digital 
image processing, light linear adjustment of brightness and 
contrast, was done using Image! 

The primary antibodies used were anti-H3K27me3 (kind 
gift of Gunnar Schotta) and anti-GFP (mixture of mouse 
monoclonal antibodies from Roche). For detection, the cells 
were stained with donkey anti-mouse Alexa 488 (Jackson 
Immuno Research), goat anti-rabbit Cy3 (Millipore), TO- 
Pro-3 (Invitrogen). For live cell imaging, MEF transfected 
cells were plated on a glass bottom p35 dish and grown 
under standard conditions. 3D stacks were acquired at 20 
min intervals overnight. Time lapse imaging was carried 



out on a Ultra VIEW VoX spinning disk confocal system 
(PerkinElmer, UK) in a closed live cell microscopy cham- 
ber (ACU control, Olympus, Japan) heated to 37°C, with 
5% CO2 and 60% air humidity control, mounted on a 
Nikon Ti microscope (Nikon, Japan), using a 60 x / 1 .45 NA 
Planapochromat oil immersion objective lens. Images were 
obtained with a cooled 14 bit EMCCD camera (C9100-50, 
CamLink). Maximum intensity Z-projections were assem- 
bled into videos and annotated using ImageJ (http://rsb. 
info.nih.gov/ij/). For cyclobutane-pyrimidine dimer (CPD) 
analysis, cells were grown on glass coverslips. Fixation of 
the cells was done with 3.7% formaldehyde for 10 min at 
room temperature followed by a combined permeabiliza- 
tion and DNA denaturation step with ice-cold 0.7% Triton 
X- 100/0.4 N HCL/PBS for 10 min. For blocking step, the 
cells were incubated with 4% BSA for 30 min. Immunos- 
taining for CPD was performed by incubating the cells with 
monoclonal anti-CPD (Kamiya Biomedical Company, Cat 
no.: MC-062) at a dilution of 1:200 for 1 h at room temper- 
ature. For detection, cells were incubated with polyclonal 
anti-mouse IgG Cy3 (Jackson Laboratories, Cat no.: 715— 
165-151) at a dilution of 1:330 for 1 h at room tempera- 
ture. The DNA content was stained with DAPI and finally 
mounted in Mowiol 4-88 (Sigma- Aldrich Chemie). 
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Figure 3. Chromatin containing the GFP-tagged histone variants shows different protein composition. (A) Coomassie blue staining of SDS-PAGE gel of 
chromatin fibers purified by GFP-nanotrap beads. Chromatin input and purified fractions are shown. The bands corresponding to the GFP-tagged histones 
are indicated: (*) GFP-H2A.Bbd, (**) GFP-H2A and (***) GFP-macroH2A.1.2. (B) Proteomic analysis of GFP-H2A.Bbd-, GFP-macroH2A.1.2- and 
GFP-H2A-containing chromatin. Specific proteins enriched in the GFP-H2A.Bbd pull-down over the proteins in the GFP-macroH2A.1.2 (x-axis) and 
GFP-H2A (y-axis) pull-downs are displayed in a scatter plot according to their log2 enrichment values. 



In situ replication labeling 

To visualize newly synthesized DNA, cells grown on cover- 
slips were pulse labeled with either 10 |xM EdU (Figure 4) or 
100 |xM 5-bromo-2 / -deoxyuridine (BrdU, Sigma- Aldrich) 
(Figure 5) and either fixed immediately or chased for the 
indicated time in pre-warmed medium supplemented with 
200 |xM thymidine. Cells were fixed with 3.7% formaldehyde 
for 10 min and permeabilized with 0.5% Triton X-100 for 
20 min. BrdU detection was performed with rat-anti-BrdU 
(Gentaur, Cat no.: OBT0030CX, dilution 1 : 100) in the pres- 



ence of DNase I as described (50). Click reaction was per- 
formed using the click-it kit (Invitrogen). Alexa Fluor 647 
azide was added in a final concentration of (2 (xM). The 
click reaction was performed for 30 min at room tempera- 
ture. The cells were then stained with an anti-GFP antibody 
as described above. 

Nuclear extract for western blot analysis 

GFP-histone variant transfected MEFs and HL cell lines 
were incubated with 0.1% Triton in PBS supplemented with 
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Figure 4. H2A.Bbd localizes to DNA synthesis sites and is replaced within 1 h. (A) Single confocal slice of female MEFs expressing the GFP-tagged 
histone variants pulsed labeled with EdU. (B) Single confocal slice of female MEFs expressing the GFP-tagged histone variants pulse labeled with EdU 
and chased for 1 h before fixation. DAPI: DNA (blue); GFP: GFP-tagged histone variant (green); EdU: replication foci (red). Scale bar: 5 |xm. The channel 
intensity recorded along the white line is plotted on a diagram. Blue line: DNA; green line: GFP-tagged histone variants; red line: replication foci. 
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Figure 5. H2A.Bbd follows PCNA at replication foci. Selected frames of 
time-lapse imaging (see Supplementary movie) of MEF cells stably ex- 
pressing GFP-H2A.Bbd (green) transfected with mRFP-PCNA (red). Z- 
stacks were collected every 20 min over night using a spinning disk micro- 
scope equipped with temperature, humidity and CO2 control. Shown are 
maximum intensity projections. Scale bars: 5 |jim. 

proteinase inhibitor for 10 min at 4°C. After centrifugation 
at 3000 rpm the nuclear pellet was resuspended in 1 x PBS 
supplemented with proteinase inhibitors and the DNA was 
digested with Benzonase (VWR) for 15 min at 37°C. 2x 
Laemmli buffer was added to the suspension of nuclei and 
boiled for 15 min at 95° C. 

Western blot antibodies 

The anti-GFP antibody (a mixture of two mouse anti-GFP 
monoclonal antibodies (Roche)) was used at a dilution of 
1:1000, the anti-H3 antibody (Abeam abl719) was used 
at a concentration of 1:5000 and the H2A.Bbd antibody 
(Millipore polyclonal) at a concentration of 1:1000. HRP- 
coupled secondary antibodies were used at a concentration 
of 1:5000 and detection was performed using ECL chemi- 
luminescent reagents (GE Healthcare Life Sciences). 

Microirradiation experiments 

Microirradiation experiments were carried out in Nikon Ti 
spinning disk microscope (see live cell imaging) with a 405 
nm diode laser set to 100% transmission. Approximately 1 
|xm spots were selected inside the nucleus and microirradi- 
ated with 1 mJ. Pre- and post-bleach images were acquired 
as indicated. 

UV-C irradiation experiments 

For UV-C irradiation experiments, the cells were irradiated 
with a UV-C source containing 5.8 W short wave lamps 
(Amersham Biosciences) at a wavelength of 254 nm. For 



damage induction, the media were removed from the cells 
and were globally irradiated for 10 s with a final dose of 1 
mJ/cm 2 . Unirradiated cells were used as control. 

CPD quantification 

The levels of DNA damage were calculated from confocal 
3D images of CPD stainings and DNA counterstained with 
DAPI using the image analysis platform Priithon (http:// 
code.google.eom/p/priithon/). The DAPI images were used 
to create 3D nuclear masks by applying a median filter 
(sigma = 500 nm) and segmenting the filtered images using 
automatic thresholding. The integrated CPD signal within 
the nuclear mask was divided by the mask volume to obtain 
the mean CPD intensity in each nucleus. 

Quantitative reverse transcription-polymerase chain reaction 

RNA from HL cell lines and PBMCs was isolated using 
Trizol reagent (Invitrogen, Karlsruhe, Germany) following 
the manufacturer's protocol. After reverse transcription of 2 
fxg of RNA, quantitative reverse transcription-polymerase 
chain reaction (qRT-PCR) was performed essentially as de- 
scribed (57). The following primer combinations have been 
used: actin beta (ACTB): 5'-ggc ate gtg atg gac tec g-3' and 
5'-gct gga agg tgg aca gcg a-3'; H2A.Bbd: 5'-tcg ttt tea gtg 
age cag gt-3' and 5'-cag aat taa tga agg ccc aag-3'. Quantifi- 
cation of gene expression was performed using the 2 _AACt 
method [H], and for comparative analysis the median Dct 
value from PBMC (21 independent donors) was set as one. 

RESULTS 

To investigate the composition of H2A.Bbd-containing 
chromatin we generated female MEF clonal cell lines ex- 
pressing GFP-tagged H2A.Bbd, macroH2A.1.2 and -H2A. 
These cell lines express the corresponding histones at ap- 
proximately similar levels (Figure 1A). After several weeks 
in culture, we noticed that whereas the levels of GFP- 
H2A.Bbd and GFP-H2A remained the same, the amount 
of GFP-macroH2A.1.2 dropped, as seen by immunoblot 
of nuclear extracts (Figure IB). This phenomenon was ob- 
served in several attempts and may be caused by a negative 
effect of high levels of macroH2A.1.2 on the cell cycle (18) 
and a subsequent negative selection of the highly expressing 
cell clones. In agreement with previously published results, 
we detect GFP-macroH2A.1.2 co-localizing with the inac- 
tive X-chromosome labeled with an antibody recognizing 
H3K27me3 (22) whereas GFP-H2A.Bbd is largely excluded 
from the Barr Body (31) (Figure 1C). The canonical GFP- 
H2A in contrast does not show enrichment on a particular 
chromosome or nuclear domain (data not shown). 

Upon establishment of stable lines expressing the differ- 
ent GFP-tagged histones, we isolated chromatin fragments 
from MEFs using micrococcal nuclease (MNase) digestion. 
As chromatin is organized in large domains in which not 
all nucleosomes contain the same histone isoform (36), we 
pooled chromatin fragments of different lengths that were 
released from the nuclei using a combination of different 
extraction protocols (Figure 2A). The tagged histone vari- 
ants (Figure 2B) were released with similar efficiency from 
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the nuclei, suggesting that the physicochemical properties of 
the nucleosomal arrays carrying specific variants in MEFs 
are roughly the same. MNase treatment released more than 
80% of the tagged chromatin, which is comparable to the 
amount isolated by Tolstorukov et al (29) and significantly 
higher than the levels released in experiments performed by 
Ishibashi et al (34). This difference may be due to experi- 
mental differences or differences in cell culture conditions 
as we also observe a small fraction of GFP-H2A.Bbd in the 
nuclear matrix fraction. However, our isolation procedure 
was selective for chromatin-bound histones, as we could not 
release histones from nuclei without MNase digestion (Fig- 
ure 2B). 

The nucleosomal arrays containing GFP-H2A.Bbd, 
GFP-macroH2A.1.2 or GFP-H2A were purified by GFP- 
affinity chromatography (Figure 3 A). The analysis of the 
protein composition of these fibres showed that, not sur- 
prisingly, the oligonucleosomal arrays contain not only the 
tagged versions of H2A used for purification but also all 
four canonical histones. This is in good agreement with our 
previous findings that chromatin domains are not homoge- 
nous with regard to their histone isoform/ variant composi- 
tion (36). 

The non-histone proteins that were isolated together with 
the tagged short nucleosomal arrays were identified and 
quantified by separation of the proteins via SDS-PAGE, 
and a subsequent identification by liquid chromatography- 
mass spectrometry (LC-MS). Proteins were identified 
from these runs using the Andromeda algorithm of the 
MaxQuant protein analysis package and quantified using 
their iBAQ values (51,58). We identified a total number of 
1428 proteins in three independent biological replicates for 
each variant and control. More than 50% of the proteins are 
identified in at least two independent biological replicates. 
The overlap varies between samples: 53.4% for H2A.Bbd, 
51.6% for H2A and 53.1% for -macroH2A. 1 .2 (Supplemen- 
tary Figure SI). A protein was considered to be specifically 
chromatin-bound when its average log2 transformed iBAQ 
value had a difference of at least 3.5 between the GFP-H2A 
variant sample and the background control and it was de- 
tected in at least two samples. The application of these filters 
resulted in the identification of 607 bona fide chromatin pro- 
teins. For every chromatin-bound protein, we determined 
its variant specificity value as the log2 value of the ratio of 
the corresponding enrichment factors (Figure 3B). As ex- 
pected based on the different chromosomal distributions 
((22,31) and Figure 1C), the GFP-H2A.Bbd- and GFP- 
macroH2A. 1 .2-containing short nucleosomal arrays were 
the most divergent (Figure 3B). For the functional anal- 
ysis, we focused our further analysis towards the proteins 
enriched in the GFP-H2A.Bbd-containing chromatin. Very 
little is known about the effect of H2A.Bbd on chromatin 
structure and its molecular function. We therefore won- 
dered whether a global analysis of the H2A.Bbd bound pro- 
teome would provide a deeper insight. To this end, we per- 
formed a gene ontology analysis for factors enriched in the 
GFP-H2A.Bbd compared to the GFP-macroH2A. 1 .2 and 
the GFP-H2A samples (upper right panel, Figure 3B). Con- 
sistent with the recent finding that H2A.Bbd might play a 
role in RNA splicing (29,30), we observe a strong overrepre- 
sentation of gene ontology categories related to RNA mat- 



uration (Tables 1 and 2). In addition, when we compared 
the H2A.Bbd- with macroH2A.l. 2-containing chromatin, 
we found a strong enrichment of factors involved in DNA 
replication such as PCNA, the MCM proteins, the FACT 
complex and the CAF1 complex (Tables 1 and 3). These 
proteins are also strongly enriched when comparing GFP- 
H2A.Bbd chromatin with chromatin isolated from cells ex- 
pressing the GFP-tagged H2A isoform (Tables 2 and 3), 
suggesting that they are specific for GFP-H2A.Bbd rather 
than excluded from chromatin that contains high levels of 
GFP-macroH2A.1.2. 

When characterizing the nuclear distribution of GFP- 
H2A.Bbd, we noticed a subpopulation of cells that had a 
particular GFP-H2A.Bbd localization pattern. In contrast 
with the hypothesis that H2A.Bbd is preferentially found 
in transcriptionally active chromatin, these cells showed a 
strong association with DAPI dense regions within the nu- 
cleus, which are usually considered to be transcriptionally 
inert (Supplemental Figure S2). Based on our observation 
that GFP-H2A.Bbd is associated with the replication ma- 
chinery, we wondered whether the cells that showed this par- 
ticular localization pattern are in fact engaged in replicating 
heterochromatic DNA. 

We therefore labeled asynchronously growing MEFs ec- 
topically expressing GFP-H2A.Bbd with EdU to identify 
cells that are in S-phase. Consistent with our hypothesis we 
observed a striking co-localization of the EdU signal with 
GFP-H2A.Bbd in all S-phase cells (Figure 4A and Supple- 
mentary Figure S3). However, after a 1 h chase following the 
30 min EdU labelling, the co-localization is lost suggesting 
a transient interaction (Figure 4B). To get a better under- 
standing of the dynamic localization of H2A.Bbd to sites of 
DNA synthesis, we performed live cell imaging of cell lines 
expressing GFP-tagged H2A.Bbd (Figure 5 and Supple- 
mentary movie). High-resolution microscopy revealed that 
H2A.Bbd is incorporated slightly after PCNA is recruited 
to sites of DNA replication, which confirms our finding that 
H2A.Bbd interacts with PCNA-containing chromatin (Fig- 
ure 5). Using this technique we could also determine a res- 
idence time of H2A.Bbd at replication foci of less than 10 
min at early replicating foci and 30 min at sites of late repli- 
cation (Supplementary Figure S4). As we see a close asso- 
ciation with PCNA, we wondered whether GFP-H2A.Bbd 
is also incorporated at sites of DNA damage. For this we 
microirradiated MEFs expressing GFP-tagged H2A.Bbd at 
450 nm for 1200 ms. After irradiation, PCNA quickly local- 
izes to sites of DNA damage followed by GFP-H2A.Bbd, 
which is similar to their timing of localization to DNA repli- 
cation sites. However, in contrast to replication sites, GFP- 
H2A.Bbd can be detected at repair foci for at least 2 h (Fig- 
ure 6). 

In order to investigate whether the expression of 
H2A.Bbd has an effect on replication and DNA damage, 
we compared the S-phase length of cells expressing GFP- 
H2A.Bbd with the ones expressing the canonical GFP- 
H2A. In agreement with the hypothesis that H2A.Bbd fa- 
cilitates chromatin remodeling or assembly during DNA 
replication, we detect a 20% shorter S-phase in the GFP- 
H2A.Bbd expressing MEFs compared to the GFP-H2A 
expressing cells (Figure 7A). Moreover, the expression of 
H2A.Bbd also increases the sensitivity to UV based DNA 
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Table 1. DNA replication factors are significantly over-represented in the GFP-H2A.Bbd compared to the GFP-macroH2A.1.2 pull-down 



Category 


Term 


Count 


0/ 
70 


77-value 


Bonferroni 


Benjamini 


1 UK 


GOTERM BP FAT 


GO:0006396-RNA processing 


100 


17.5 


3.32E-59 


3.71 E-56 


3.71 E-56 


5.33E-56 


GOTERM BP FAT 


GO:0006412~translation 


75 


13.2 


1.35E-44 


1.51E-41 


7.55E-42 


2.17E-41 


GOTERM BP FAT 


GO0006397~mRNA nrocessine 


65 


11.4 


5.53E-40 


6.17E-37 


2 06E-37 


8 86E-37 


GOTERM BP FAT 


GO:0042254~ribosome biogenesis 


45 


7.9 


1.22E-37 


1.36E-34 


3.40E-35 


1.95E-34 


GOTERM BP FAT 


GO:0008380-RNA splicing 


56 


9.8 


3.24E-37 


3.62E-34 


7.23E-35 


5.19E-34 


GOTERM BP FAT 


GO:0016071~mRNA metabolic 


66 


11.6 


5.59E-37 


6.24E-34 


1.04E-34 


8.96E-34 




process 














GOTERM BP FAT 


GO:002261 3~ribonucleoprotein 


47 


8.2 


1.10E-35 


1.23E-32 


1.76E-33 


1.76E-32 




complex biogenesis 














GOTERM BP FAT 


GO:0006364~rRNA processing 


29 


5.1 


9.06E-24 


1.01E-20 


1.26E-21 


1.45E-20 


GOTERM BP FAT 


GO:0016072~rRNA metabolic 


29 


5.1 


1.40E-23 


1.57E-20 


1.74E-21 


2.25E-20 




process 














GOTERM BP FAT 


GO 1 000632 5 ^chromatin organization 


47 


8.2 


5.21 E-l 9 


5.82E-16 


5.82E-17 


8.35E-16 


GOTERM BP FAT 


GO:005 1 276~chromosome 


53 


9.3 


6.67E-19 


7.44E-16 


6.77E-17 


1.07E-15 




organization 














GOTERM BP FAT 


GO:0034470~ncRNA processing 


31 


5.4 


7.42E-16 


8.67E-13 


7.23E-14 


1.24E-12 


GOTERM BP FAT 


GO*0034660~ncRNA metabolic 


31 


5.4 


6.50E-13 


7.25E-10 


5.58E-11 


1.04E-09 




process 














GOTERM BP FAT 


GO0006260-DNA renlication 


27 


4.7 


8.31 E-l 3 


9.28E-10 


6.63E-11 


1.33E-09 


GOTERM BP FAT 


GO:0016568~chromatin modification 


33 


5.8 


1.38E-12 


1.54E-09 


1.03E-10 


2.22E-09 


GOTERM BP FAT 


GO:0006259-DNA metabolic process 


44 


7.7 


3.61 E-12 


4.03E-09 


2.52E-10 


5.78E-09 


GOTERM BP FAT 


GO:0016569~covalent chromatin 


20 


3.5 


7.62E-11 


8.50E-08 


5.00E-09 


1.22E-07 




modification 














GOTERM BP FAT 


GO:0016570~histone modification 


19 


3.3 


2.95E-10 


3.30E-07 


1.83E-08 


4.73E-07 


GOTERM BP FAT 


GO:0006333~chromatin assembly or 


20 


3.5 


7.51 E-l 0 


8.38E-07 


4.41 E-08 


1.20E-06 




disassembly 














GOTERM BP FAT 


GO:0034622~cellular 


25 


4.4 


5.55E-08 


6.20E-05 


3.10E-06 


8.89E-05 




macromolecular complex assembly 














GOTERM BP FAT 


GO:0016573~histone acetylation 


11 


1.9 


6.14E-08 


6.86E-05 


3.27E-06 


9.84E-05 


GOTERM BP FAT 


GO:0006473~protein amino acid 


11 


1.9 


8.28E-08 


9.24E-05 


4.20E-06 


1.33E-04 




acetylation 














GOTERM _BP_FAT 


GO:003462 1 -cellular macromolecular 


25 


4.4 


5.41 E-07 


6.04E-04 


2.63E-05 


8.67E-04 




complex subunit organization 














GOTERM BP FAT 


GO:0006350~transcription 


91 


16.0 


6.46E-07 


7.21 E-04 


3.00E-05 


















0.0010346 


GOTERM BP FAT 


GO:0043543~protein amino acid 


11 


1.9 


1.53E-06 




6.84E-05 


0.0024541 




acylation 








0.00170809 






GOTERM BP FAT 


GO:0043967-histone H4 acetylation 


8 


1.4 


1.95E-06 


0.0021771 


8.38E-05 





0.00312867 



Top enrichment terms for the factors enriched in the H2A.Bbd/macroH2A.1.2 sample. The terms and their related /^-values for Biological Processes (BP) 
are displayed. Category: original database/resource where the terms orient; Term: enriched term associated with the gene list; Count: gene in the list 
involved in the term; Percentage: involved genes/total genes; />-value: modified Fisher exact /?-value, EASE score. 



damage as measured by the increased formation of CPDs 
(Figure 7B). The shortening of S-phase and increased sen- 
sitivity toward DNA damage are hallmarks of cancer cells. 
Indeed, the analysis of Gene Expression Omnibus microar- 
ray data sets from HL cell lines and primary HL samples in 
comparison with normal tissue revealed increased levels of 
H2A.Bbd RNA in a large number of HL lines ((17); Fig- 
ure 7C). Consistent with this finding, chromatin isolated 
from HL cells contains much higher levels of H2A.Bbd 
than control cells (Figure 7D) where H2A.Bbd is not de- 
tectable. Interestingly, the total levels of H2A.Bbd vary sub- 
stantially in different HL lines with the ones containing a lot 
of H2A.Bbd having a shorter doubling time than the ones 
with lower expression levels ((32); Figure 7E), suggesting 
that H2A.Bbd overexpression might have a similar effect in 
MEFs and tumor cells. 

DISCUSSION 

So far most studies addressing the function of the histone 
variant H2A.Bbd focused either on the structural character- 



istics of nucleosomes containing these variants (33) or on 
the DNA sequences that are associated with it (27,29,30). 
However, the physicochemical properties of chromatin and 
its function are mainly defined by associated proteins. As 
the current knowledge about the proteome bound to the 
H2A.Bbd-containing chromatin is limited, we performed a 
detailed proteomic analysis on chromatin isolated from cells 
expressing a GFP-tagged version of H2A.Bbd and used 
the purification of nucleosomal arrays containing GFP- 
macroH2A. 1 .2 and GFP-H2A as highly specific filters for a 
differential proteomics approach. As the cells we use to an- 
alyze H2A.Bbd bound chromatin do not show significant 
expression levels of this variant, the system most likely re- 
flects the situation in tumor cells that have a high expression 
level for H2A.Bbd. 

In agreement with earlier findings investigating chro- 
matin fibres containing H3 variants H3.1 and H3.3 (36), 
we observed that arrays of nucleosomes containing either 
one of the two H2A variants also contain a sizeable amount 
of canonical histones. Despite this apparent heterogeneity 



Nucleic Acids Research, 2014, Vol 42, No. 10 6415 



Table 2. Replication factors are significantly over-represented in the GFP-H2A.Bbd compared to the GFP-H2A pull-down 



Category 


Term 


Count 


% 


Rvalue 


Bonferroni 


Benjamini 


FDR 


GOTERM BP FAT 


GO:0006396-RNA processing 


90 


18.6 


1.61E-55 


1.66E-52 


1.66E-52 


2.56E-52 


GOTERM BP FAT 


GO:0006412~translation 


73 


15.1 


7.94E-48 


8.15E-45 


4.07E-45 


1.26E-44 


GOTERM BP FAT 


GO:0006397~mRNA 


60 


12.4 


5.30E-39 


5.43E-36 


1.81E-36 


8.40E-36 




processing 














GOTERM BP FAT 


GO:0008380-RNA splicing 


53 


11.0 


1.18E-37 


1.21E-34 


3.02E-35 


1.87E-34 


GOTERM BP FAT 


GO:0016071~mRNA metabolic 


61 


12.6 


2.38E-36 


2.44E-33 


4.88E-34 


3.77E-33 




process 














GOTERM BP FAT 


GO:0042254~ribosome 


40 


8.3 


6.15E-34 


6.31 E-31 


1.05E-31 


9.75E-31 




biogenesis 














GOTERM BP FAT 


GO :00226 1 3 ~ribonucleoprotein 


41 


8.5 


2.40E-31 


2.46E-28 


3.52E-29 


3.80E-28 




complex biogenesis 














GOTERM BP FAT 


GO:0006364~rRNA processing 


24 


5.0 


3.70E-19 


3.80E-16 


4.75E-17 


5.87E-16 


GOTERM BP FAT 


GO:0016072~rRNA metabolic 


24 


5.0 


5.21 E-l 9 


5.35E-16 


5.94E-17 


8.26E-16 




process 














GOTERM BP FAT 


GO:0034470~ncRNA 


26 


5.4 


3.25E-13 


3.34E-10 


3.34E-11 


5.16E-10 




processing 














GOTERM BP FAT 


GO:0034660~ncRNA 


26 


5.4 


8.46E-11 


8 68E-08 


7.89E-09 


1.34E-07 




metabolic process 














GOTERM BP FAT 


GO:0006325~chromatin 


31 


6.4 


7.43E-10 


7.62E-07 


6.35E-08 


1.18E-06 




organization 














GOTERM BP FAT 


GO:005 1 276~chromosome 


34 


7.0 


5.25E-09 


5.39E-06 


4.14E-07 


8.32E-06 




organization 














GOTERM BP FAT 


GO:0006260-DNA replication 


20 


4.1 


1.39E-08 


1.42E-05 


1.02E-06 


2.20E-05 


GOTERM BP FAT 


GO:0006259~DNA metabolic 


31 


6.4 


5.34E-07 


5.48E-04 


3.66E-05 


8.47E-04 




process 














GOTERM BP FAT 


GO:0016568~chromatin 


22 


4.5 


8.78E-07 


9.01 E-04 


5.63E-05 






modification 












0.00139236 


GOTERM BP FAT 


GO:0034622~cellular 


21 


4.3 


9.26E-07 


9.50E-04 


5.59E-05 





macromolecular complex 0.00146848 
assembly 



Top enrichment terms for the factors enriched in the H2A.Bbd/H2A sample. The terms and their related /7-vames for Biological Processes (BP) are 
displayed. Category: original database/resource where the terms orient; Term: enriched term associated with the gene list; Count: gene in the list involved 
in the term; Percentage: involved genes/total genes; /7-value: modified Fisher exact /7-value, EASE score. 



Table 3. Replication associated factors are enriched in chromatin containing GFP-H2A.Bbd 



Iog2 enrichment 

Iog2 enrichment GFP-H2A.Bbd/GFP- 
Protein names GFPH2A.Bbd/GFP-H2A /?-value macroH2A.1.2 /?-value 



Pcna 


9.26 


0.00555361 


9.03 


0.0041637 


Ginsl 


7.42 


0.01801521 


5.80 


0.0312951 


Mcm3 


5.44 


0.2943226 


4.84 


0.2709437 


Cdc21;Mcm4 


-1.04 


0.4773969 


0.09 


0.9612214 


Cdc46;Mcm5 


0.25 


0.9734912 


6.21 


0.1672788 


Cdc47;Mcm7 


1.17 


0.7511821 


2.63 


0.2726894 


Mcm6 


0.86 


0.9215584 


3.21 


0.4707804 


Ssrpl 


4.31 


0.1783074 


10.13 


0.0091877 


Fact 140 


5.35 


0.1497874 


9.86 


0.0166473 


Chaflb 


7.02 


0.09077201 


5.98 


0.1004223 


Rbbp4 


2.56 


0.6075331 


3.46 


0.3578618 


Rfc4 


8.25 


0.1551118 


9.51 


0.0762587 


Rfc2 


6.89 


0.1028537 


7.39 


0.0611908 


Rfc3 


1.49 


0.9385401 


2.87 


0.8451195 



Enriched replication factors in the GFP-H2A.Bbd-containing chromatin are displayed according to their enrichment factors and p- values calculated in 
the GFP-H2A.Bbd/GFP-macroH2A.1.2 and in the GFP-H2A.Bbd/GFP-H2A samples. 



of the purified fibres, these histone variants have a major 
impact on the proteomic composition of chromatin. Note- 
worthy, the number of proteins detected in GFP-H2A.Bbd 
chromatin is higher compared to chromatin isolated from 
GFP-H2A or GFP-macroH2A 1.2. This observation could 
be due to an intrinsic feature of H2A.Bbd containing chro- 
matin or a result of slight expression differences of the his- 
tone variants within the cells. In order not to bias our quan- 



titative analysis, we normalized the resulting quantitative 
information of all variant specific chromatin proteins to the 
total amount of isolated chromatin using the histone H4 
iBAQ values. 

We observed a striking enrichment of factors involved in 
DNA replication within the GFP-H2A.Bbd-containing mi- 
crosomal arrays in an asynchronous population of cells. 
Although most proteins reported to reside in the mam- 
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Figure 6. H2A.Bbd accumulates at DNA damage sites following PCNA. Time-lapse imaging of microirradiated MEF cells stably expressing GFP- 
H2A.Bbd (green) and transfected with mRFP-PCNA (red) at 405 nm for 1200 ms. Images were collected before microirradiation and 5 s, 10/50 and 
110 min after microirradiation. Scale bars: 5 |jim. 



malian replisome are enriched in the GFP-H2A.Bbd asso- 
ciated chromatin, not all enrichment values are significantly 
different among the variant containing chromatin when an- 
alyzed by rigorous statistical testing. This may well be due 
to the fact that only a fraction of cells within the popula- 
tion is actively replicating and incorporating H2A.Bbd at 
sites of DNA synthesis. 

Nevertheless, the presence of factors associated with 
DNA replication in the chromatin containing GFP- 
H2A.Bbd prompted us to analyze its involvement in this 
process using alternative methods. Indeed, when we ana- 
lyzed DNA replication by immunohistochemistry and live 
cell imaging, we observed a strong co-localization pattern 
of GFP-H2A.Bbd with sites of DNA replication. Interest- 
ingly, a previous report had reported a dotted staining pat- 
tern in a small fraction of H2A.Bbd transfected MEFs (34), 
which in light of our findings most likely reflects replicat- 
ing chromatin. Histone variants have been shown to be re- 
cruited to DNA repair foci (37), however, they have not 
been found co-localizing with replicating DNA. H2A.Bbd 
is a striking exception to this rule. Interestingly, this enrich- 
ment is only transient as the co-localization of H2A.Bbd 
with newly synthesized DNA is only detected for less than 
60 min, which may be the reason why it has gone unnoticed 
so far. This apparent instability of the newly formed nucle- 
osome is in agreement with the observed low intrinsic sta- 
bility of the H2A.Bbd-containing particles and its inability 
to bind the linker histone HI (39,41). 

H2A.Bbd is incorporated at both early and late replicat- 
ing chromatin during S-phase, which points to a general role 
of the DNA synthesis machinery in the incorporation of this 
histone variant. Alternatively, the localization of H2A.Bbd 
to sites of DNA synthesis might reflect its ability to bind to 



naked DNA at replication sites. However, we do not think 
this is the sole mechanism of incorporation as we do not ob- 
serve it for the canonical H2A variant that is expressed to 
comparable levels and has a similar affinity for naked DNA. 

The differential residence time of H2A.Bbd in early or 
late replicating foci may be related to the different efficiency 
of replication fork progression in eu- and heterochromatin 
or to a different binding dynamics of an H2A.Bbd inter- 
acting factor. What is intriguing, however, is the substantial 
increase in residence time on DNA repair foci, which may 
be due to an extended presence of PCNA at those regions. 

Like most other studies performed before, we do not in- 
vestigate H2A.Bbd containing chromatin isolated from cells 
that naturally express this variant. It is therefore difficult to 
make a conclusion about its physiological function. How- 
ever, it is tempting to speculate that the incorporation of 
H2A.Bbd at sites of DNA synthesis could physiologically 
contribute to the maturation of sperm and the replace- 
ment of the canonical histone with protamines. The puta- 
tive mouse homologue H2A.Lapl is expressed in all stages 
of spermatogenesis with the highest level between the late 
pachytene and round spermatid stage when it is exported to 
the cytosol and another H2A.Bbd-like histone becomes in- 
corporated into the chromatin (17-20,35). The presence of 
H2A.Bbd-like histones in mouse coincides with episodes of 
increased DNA breakage in sperm (17,18,21,35,59). More- 
over, the DNA repair machinery might facilitate the incor- 
poration of H2A.Bbd, thereby facilitating DNA recombi- 
nation and protamine exchange. As spermatocytes express- 
ing H2A.Bbd are postmitotic, it is unlikely that H2A.Bbd 
plays a role in DNA replication. 

As pointed out above, our system more closely resembles 
the elevated expression of H2A.Bbd in HL cells (Figure 7C 
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Figure 7. MEF cells and HL lines expressing H2A.Bbd replicate faster and are more sensitive to DNA damage. (A) Doubling time and S-phase length in 
MEF cells stably expressing GFP-H2A.Bbd and GFP-H2A. (B) Mean intensity of nuclear CPD signal 10 and 30 min post 1 mJ UV-C irradiation. Error 
bars represent the 95% confidence interval. N > 30 (C) qRT-PCR expression analysis of H2A.Bbd in HL cells and control cells. cDNA from HL cell lines 
(two independent experiments) and normal PBMC (21 donors) was used as template for quantitative PCR. For calculation of relative expression values, 
actin beta was used as housekeeping control and the median expression in normal PBMC was set as 1. (D) Western blot analysis of H2A.Bbd in Hodgkin's 
lymphoma cells (HL) and control cell lines. Ponceau staining of histone H3 is used as loading control. (E) Dot plot showing the correlation in Hodgkin's 
lymphoma cell lines between the doubling time derived from (32) (y-axis) and the expression levels of H2A.Bbd (x-axis) measured by qRT-PCR. Every dot 
represents a specific Hodgkin's lymphoma cell line. 
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and D). Such an overexpression might result in a faster S- 
phase progression similar to what we observe in MEFs ex- 
pressing GFP-H2A.Bbd (Figure 7). Indeed, HL cell lines 
with lower expression of H2A.Bbd (cell lines L-540 and 
HDLM-2) have longer doubling times than cell lines with 
higher expression of H2A.Bbd (22,23,32). HL high sensi- 
tivity to radio- and chemotherapy is known for a long time 
and has led to the development of highly effective treatment 
regimens (reviewed in (60)). H2A.Bbd is specifically overex- 
pressed in HL cell lines and in HL cells in vivo (17) and might 
be involved in the high success of radio- and chemotherapy 
for treatment of these patients. 

CONCLUSION 

Our data support a functional compartmentalization of 
chromatin by histone variants and their associated factors. 
The finding that H2A.Bbd transiently localizes to sites of 
DNA synthesis in cells thereby increasing S-phase progres- 
sion and sensitizing the cells for DNA damaging UV irra- 
diation opens up new avenues for the investigation of the 
function of this exceptional histone variant in healthy and 
diseased tissue. 
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